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Abstract 

A new approach to the theory of anisotropic exciton based on Fock transformation, i.e., on a 
stereographic projection of the momentum to the unit 4-dimensional (4D) sphere, is developed. 
Hyperspherical functions are used as a basis of the perturbation theory. The binding energies, 
wave functions and oscillator strengths of elongated as well as flattened excitons are obtained 
numerically. It is shown that with an increase of the anisotropy degree the oscillator strengths 
are markedly redistributed between optically active and formerly inactive states, making the latter 
optically active. An approximate analytical solution of the anisotropic exciton problem taking 
into account the angular momentum conserving terms is obtained. This solution gives the binding 
energies of moderately anisotropic exciton with a good accuracy and provides a useful qualitative 
description of the energy level evolution. 
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I. INTRODUCTION 



The interest to the anisotropic exciton problem^ has been revived with the progress 
in the physics of semiconductor heterostructures. In semiconductor superlattices the mini- 
band formation causes a strong mass anisotropy.- In fact, the localization of carriers inside 
quantum wells and their tunneling trough barriers can be described in terms of anisotropic 
medium approximation as the effect of mass renormalization. The dielectric constant be- 
comes anisotropic also if the superlattice constituent layers have different dielectric suscep- 
tibilities. Recently such a formalism has been used in the theory of excitons in short-period 
superlattices (see, e.g., Refs. |5|j6|). 

The main complication of the uniaxial anisotropic exciton problem is that the Coulomb 
potential symmetry is broken (the spherical symmetry as well as the "hidden" one, the in- 
trinsic property of the hydrogen-like system) so that only the angular momentum projection 
and parity conserve. As a consequence, the solution of the Schrodinger equation is no more 
factorized into radial and angular parts and cannot be represented as a finite combination 
of standard special functions. 

The anisotropic exciton problem was first studied by Kohn and Luttinger- (for donor 
states in silicon and germanium) by means of the variational approach with allowance for 
a group symmetry of the particular materials. Further theoretical studies^— were focused 
on perturbative solutions of the anisotropic exciton problem. For slightly anisotropic sys- 
tem Hopfield and Thomas 7 found the first-order solution, treating the anisotropy of the 
kinetic energy as a perturbation™ linear in the anisotropy parameter. The effects in a weak 
magnetic field also have been taken into account in this approximation. For a moderate 
exciton anisotropy Wheeler and Dimmock- used an expansion of the anisotropic potential 
over its asymmetric part z 2 /r 2 up to the second order in the anisotropy parameter terms, 
thus calculating in part the second-order perturbation solution. This partial diagonalization 
was completed by Deverin,- who considered the diagonal elements of the exact anisotropic 
kinetic energy (for nondegenerate levels) as well as the transcendental solution of a secular 
problem for degenerate levels. The full expansion of the anisotropic potential was con- 
sidered by SegaU^, where only the spherically symmetric part of the full expansion was 
taken into account. Finally, Faulkner 1 ^ performed calculations of donor energy levels by 
means of Rayleigh-Ritz perturbation method containing numerous (depending on hydrogen 
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quantum numbers) variational parameters. Being included in the radial part of hydrogen 
basis functions, these variational parameters served as scaling factors depending on the 
anisotropy degree. In the limit of an extreme anisotropy, the exciton binding energies were 
calculated^ 2 -' in adiabatic approximation. Following the method suggested by Faulkner, 
Baldereschi and Diaz' obtained similar results and attempted to calculate excitonic oscil- 
lator strengths. The same Rayleigh-Ritz method was used in Ref. [l7| for calculations of the 
energy levels of 2D anisotropic exciton. 

Recently, an elegant model of fractional- dimensional space has been developed [see 



Refs. 



19 



20l and references therein]. It allows to treat self-consistently the bound as well 



as continuum states in hydrogen problem of noninteger dimension. However, its direct 
applicability to the anisotropic exciton problem is problematic. The reason is that the 
fractional- dimensional hydrogen problem conserves the Coulomb degeneracy of levels (so 
that the binding energies depend on the principal quantum number only), whereas in reality 
the anisotropy lifts this degeneracy and restores it only in 2D and 3D cases. 

In spite of a long history of theoretical study, the investigation of the optical properties 
of the anisotropic exciton is still not complete. For example, the behavior of exciton oscilla- 
tor strengths is very important for the understanding the experimental absorption spectra. 
However, the evolution of the oscillator strengths of the anisotropic exciton with the increase 
of the anisotropy has not been investigated, for our knowledge, with two exceptions: calcula- 
tions for slightly anisotropic exciton' and simulations of optical spectra within an isotropic 
exciton model.— One should note that none of the approaches'" 2 ^ is able to describe the 
drastic changes of oscillator strengths (due to the level anticrossings') with increase of the 
anisotropy reported in our paper. 

In the present paper we develop 22 a perturbation approach to the uniaxial anisotropic 
exciton problem, based on the method of stereographic projection of the momentum space 
to the unit 4D-sphere, proposed by Fock.— We use the hyperspherical harmonics, i.e., the 
irreducible representation of rotation group 0(4) of a 4D-sphere, as a basis of Brillouin- 
Wigner perturbation method. 

This approach has a number of advantages and clarifies the physical properties of the 
anisotropic exciton. (i) It allows us to utilize the additional hidden symmetry of Coulomb 
potential for expansion of anisotropic exciton wave function. Namely, for the bound exciton 
states the irreducible representation of the full symmetry group 0(4) constitutes a complete 
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set for such expansion. This expansion depends explicitly on the exciton energy through 
scaling parameters which follow adiabatically the changes in anisotropy. These parame- 
ters, similar to those introduced in the Rayleigh-Ritz method 11 (where they were defined 
by minimizing the energy functional) are exactly determined in our method. As a result, 
the hyperspherical functions turn out to be the most effective basis for numerical calcula- 
tions, (ii) Within Fock representation, the hydrogenic spectrum with the level series limit 
transforms into an equidistant one, which provides a good convergence of our method in 
a wide region of the anisotropy parameter, (iii) The matrix elements of the perturbation 
are found as analytical elementary expressions, (iv) This analytical form of perturbation 
matrix elements allows us to construct a spherical approximation with an analytical solution 
and to summarize exactly the rest part of perturbation in the second order. This spherical 
approximation, which works well in the region of a moderate anisotropy, turns out to be 
very useful for qualitative classification of the energy levels. 

We calculate numerically the energy spectrum, excitonic wavefunctions and oscillator 
strengths for flattened as well as elongated excitons. 

The paper is organized as follows. In Sec. II the expansion is formulated on the basis 
of hyperspherical formalism and basic equations of the perturbation method are derived. 
Results and discussions are presented in Sec. III. 



II. ANISOTROPIC EXCITON IN FOCK REPRESENTATION 



A. Hyperspherical formalism 

The Hamiltonian of the uniaxial anisotropic exciton is given by 

2 t jL ± \dx 2 + dy 2 ) 2/i||fe 2 y/e\\e ± (x 2 + y 2 ) + e'iz 2 ' U 

Here /i is the reduced exciton mass, e is the semiconductor dielectric constant, and sub- 
scripts || and _L refer to the quantities along and normal to the axis of symmetry (z-axis), 
respectively. In Eq. (CQ) both the kinetic and potential energies are anisotropic. However, a 



dilatation z — > Zy/e\\/e± makes the potential energy spherically symmetric. In the effective 
atomic units 

t> * _ /^e 4 * _ h2£ o m 
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where Eq = ^ye±e\\, Eq. (GQ) takes the form 



p z + ef z -~ ) V(r) = £?V»(r). 



(3) 



Here we introduced the perturbation parameter, e = 7 — 1, connected to the anisotropy 
parameter, 
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£||/i|| 



(4) 



(0 < 7 < 1 and 1 < 7 < 00 for, respectively, flattened and elongated exciton), p and p z 
denote, respectively, the dimensionless operators of momentum and its z-projection. 

We investigate the bound states with eigenenergies E v < 0, measured in Ry*, Eq. (|2J). It 
is convenient to introduce a parameter (for each bound states v) 



P, 



v = \/—E u , 



(5) 



which will play the role of the adiabatic parameter in the perturbation theory. After the 
Fourier transform, Eq. fl3]) takes the integral form 

1 r Mp' 



2tt 2 



IP - P 



(6) 



Following Fock's paper,— we perform a stereographic projection of 3D momentum space to 
the 4D unit sphere, p/p u — > u, where the 4D vector u on the sphere is defined as 

2p u p p 2 - p\ 



u = {u,m„} 



p 2 + pi ' p 2 + pi 

p = |p|. In the hyperspherical coordinates, (a, 9, <p), the unit vector u takes the form 



(7) 



u. 



u. 



u 7 



u, 



P +PZ 



p 2 +pl 



ZPuPz 

p 2 +pI 



p 2 -pI 
p 2 +pI 



sin ot sin 9 cos ip, 



sin a sin 9 sin ip, 



sin a cos 9 , 



cos a, 



(8) 



and 



d Q = sin ada sin 



ipl 



(p 2 + pi) 



-d 3 p. 



(9) 
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Let us introduce a new wave function 

*M = {p2+ 5 /p 2 Mp), (io) 

with normalization condition 

J I^(P)| 2 ^ 3 P = J (I ~ cosa)|tt„(ii)| Vfi = 1. (11) 
Then Eq. fl6]) takes the form 

(l + e -v) m u {u) = -H^ v {u). (12) 

Here ifo is the Hamiltonian of unperturbed (hydrogen-like) problem, 

1 f W 



H V(u) = — / —^-d A Q', (13) 
2ir z / \u — u'r 



and V is the perturbation operator, 



u 2 



V = z — = (1 + cos a) cos 2 0. (14) 

1 - Mr 



If e = (or 7 = 1), Eq. ( IT2|) describes the isotropic 3D exciton. As it was shown by Fock,— 
the solutions of the integral equation 

tf<°> = A(°^o* {0) (15) 

are 



*2LM,¥>) = (-20'" J ^ + ^ oCia-i(«» ¥>), (16) 

= n = l,2,..., Z = 0,...,n-1, m = 0, ±1, . . . , ±1. (17) 

Here C™(x) are the Gegenbauer polynomials^ and Yi m (9,ip) are the conventional spherical 
harmonics. The hyperspherical functions, Eq. (fT6l) . afford the irreducible representation 
of the full symmetry group 0(4) of the hydrogen-like system.— Due to the properties of 
irreducible representations, the hyperspherical function are orthogonal and normalized as 

/ \9<SL(a,6,<p)\'*n = |(l-cosa)|vI/W m (a,0,^)| 2 ^ = l, (18) 

in accordance 2 ^ with Eq. (ITT]) . It can be shown^l that the standard hydrogen wave function 
JP) 



nlm 



'r) with a given set of quantum numbers (n,l,m) (see, e.g., in Ref.U) can be Fourier 



transformed into the hyperspherical function, Eq. f|T6|) . 



B. Formulation of Brillouin-Wigner perturbation theory 



We use the Brillouin-Wigner perturbation theory, i.e. the direct diagonalization of a 
truncated Hamiltonian matrix in order to solve the anisotropic exciton problem in the form 
of Eq. (TT2"j) . The set of the hydrogen bound states eigenfunctions is not complete and the 
scattering states also must be taken into account. However, in Fock representation we are 
able to construct a complete basis out of the set of the hydrogen bound states. As it was 
shown in Ref. |25|, the scattering states are mapped on a two-sheeted hyperboloid in a 4D 
space with Minkowski metrics, whereas the bound states are mapped into a unit sphere 
via the transformation Eq. (jSJ). Thus, the problems of the bound and scattering states 
are mapped onto different subspaces, each of them to have its own complete basis. The 
anisotropic problem is mapped into the same subspaces through the transformation Eqs. (JTJ)- 
(|T0|) for the bound states and the corresponding procedure (with positive energies) for the 
scattering states. So, being interested in bound states in the whole physical region — 1 < 
e < oo, excluding the points e = — 1 (purely 2D exciton) and e = oo (purely ID exciton), we 
can use the hyperspherical harmonics Eq. ffTB 7 ]) as a complete set of basic functions.— As it 
immediately appears from Eqs. ( II 2p and (I14p . the perturbation scheme converges for |e| < 1. 
For the opposite case of e > 1 we can reformulate the perturbation problem with the help 
of the transformation p 2 + ep 2 — (e + l)[p 2 + (l/(e + 1) — l)(p 2 + Py)}- After this, we can 
redefine the effective atomic units Eq. ([2]) and consider the operator (l/(e + 1) — l)(pl + p 2 ) 
as a perturbation, thus providing the convergence for |l/(e + 1) — 1| < 1. 

The eigenfunctions are expanded as 

m u {a,0^) = S- 1 Y,C u s ^f{a ) e^) ) 8 = {n,l,m), (19) 

s 

where normalizing constants are defined as 

Si = C n,i,m [nC u n , hm - V(n + l + l)(n-l)C: +1>l!m ] . (20) 

n,l 

Then, the Schrodinger equation takes the matrix form 

( n5 ™' + \ v ») = ( 2l ) 

s' 

where 

K = 7, = vk' (22) 
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and the perturbation matrix is 

V ss > = vW J (a, e,cp)(l + cos a) cos 2 0tf J° (a, 9, (p)d 4 tt. 
Nonvanishing matrix elements V ss ' are (see Appendix A) 

Vnn™™ = Qlm j™W + ^V( n ~ l )( n + 1 + 1)<W n> + ^y/ (n - I - l)(n + l)5 n -i n > 

with 



Qlm = 77 + 



1 - 4m 2 



2 2(2/ - l)(2/ + 3)' 



and 



with 



rrl l-2;mm [I 2 ~ ™ 2 } [g ~ I) 2 ~ ™?\ , (n - I - 1)! {W + / - 2)! , 

V„„, = \ I — 7T7^ — \ I ; — — — ; —-T^l 



(2Z + l)(2Z-3) 
-1, n'<n-2. 



(n + l)\ (n'-i + iy: 



( n 2 - Anl - I 2 + n + 31 - 2 , 



n — I 



2(n- 1) 



n = n — 1, 



2n(2Z - 1) 



x < 



n-l + 1, 



(n-l + l)(n-l + 2) 
2(n + l) : 



n = n, 



n' = n + 1, 



(23) 



(24) 



(25) 



(26) 



(27) 



0, ti' > n + 2. 
All the other matrix elements vanish. 

The perturbation method in the form of Eq. (121 j) is very convenient. First of all, the 
perturbation ep 2 is invariant with respect to rotations around the £-axis and to the trans- 
formation p — > —p. Thus, each perturbed state has a definite parity and definite magnetic 
quantum number m, and the perturbation problems Eq. (1211) can be solved separately for 
different parity and m. It implies also that the summation over s' in Eq. (|2ip and thereafter 
means that only the hydrogen states with a given parity and magnetic quantum number 
have to be taken into account. The time-conjugated states ±m are still degenerate. Sec- 
ondly, the precise form of perturbation matrix V ss > provides more rigorous selection rules. 
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Namely, only the matrix elements with [see Eqs. (I24p - (I2"T]) ] 



I, n' — n, n ± 1 
I' = < 1-2, n' < n + 1, 
1 + 2, n' > n - 1 



(28) 



v 



are nonvanishing. 

The expansion (fT9{) corresponds to the following coordinate representation of the 
anisotropic exciton wave function 



where 0„^(r) are the standard hydrogen wave functions. It follows from Eq. ( 1291) that the 
wave function of anisotropic exciton takes the form of an infinite superposition of spherical 
harmonics with radially dependent coefficients. The scaling factors p u in the wave functions 
Eq. (|29|) . which are different for different perturbed states and change adiabatically with 
e, play the role of adiabatic scaling parameters in the perturbation theory. Moreover, the 
coefficients p u are analogous to the parameters in the Rayleigh-Ritz method. In contrast 
with previous works,— where p v have been calculated variationally, in our approach they 
are strictly determined by Eq. (l22il . Finally, in spite of the energy scaling factors in the basis 
functions, the effective Hamiltonian matrix in Eq. (I2T]) is energy independent, thus allowing 
for the direct diagonalization. 

It can be seen from Eq. (121]) that in Fock representation the spectrum of the unperturbed 
problem does not have a series limit. This fact is favorable for the convergence of the 
perturbation theory. Moreover, the spectrum is equidistant with respect to the hydrogen 
principal quantum number n. The matrix Eq. ( 1241) is tridiagonal. The off-diagonal matrix 
elements Eq. fl26|) with n! ~ n are rather significant but do not exceed n/2, i.e., they are 
of the order of the magnitude of the corresponding eigenvalues of the unperturbed problem. 
The other nonzero elements, V^ n ' _2 ' mm oc n~ l , decrease rapidly for fixed n' and /, and n ^> n', 
I > 1. Thus, in numerical calculations we can take into account only the states with lower 
/, and introduce a 7-dependent upper bound for the orbital quantum number. Though 
the method provides a good convergence in a large region of 7, it does not allow to avoid 
instabilities near 7 = or 7 — > 00, where the perturbation scheme becomes unstable, and a 
strong mixing of levels occurs. 




(29) 



s 
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We would like to emphasize that presented perturbation method can be easily generalized 



25, where the method of 



for an arbitrary integer dimension D > 2 in accordance with Ref. 
stereographic projection has been expanded to higher dimensions. In particular, for D = 2 
the standard spherical harmonics Yi m (6, <p) have to be used as a basis and the operator 
(e/2)V = (e/2)(l + cos#)cos 2 <^ — as a perturbation. Here cos# = (p 2 — pl)/{p 2 + pi), 
tany? = p y /p x . 

The problem of anisotropic exciton scattering states can be approached analogously using 
hyperspherical harmonics on a two-sheeted 4D-hyperboloid as a basis for the perturbation 
problem. The eigenvalues should be defined with positive energies, instead of Eq. (|22[) . 
However, the eigenvalue problem [analogous to Eq. (I2ip ] becomes more complicated: we 
have to solve now a system of integral equations, because of dependence on continuum 
quantum numbers. 

One should note that the method of stereographic projection can be formally generalized 
for the fractional-dimensional exciton problem, the exciton binding energies coinciding with 



those obtained in Ref. 



19l . However, due to the generalized hyperspherical symmetry con- 



servation (the anisotropy parameter now appears in a role of the fractional dimensionality), 
the energy levels are Coulomb degenerate, as it was mentioned above. 



III. RESULTS AND DISCUSSIONS 



Due to the symmetry properties of uniaxial anisotropic exciton Hamiltonian, matrices 
with even and odd I as well as with different m can be diagonalized independently. In 
contrast to the variational technique which provides only the upper bound of the binding 
energies, the Brillouin-Wigner perturbation method allows us to reach necessary precision 
by choosing a sufficiently large matrix to be diagonalized. We perform our calculation with 
a relative energy precision of 10~ 4 . In order to provide this precision in the calculation of the 
ground state energy for 0.6 < 7 1 / 3 < 2, hydrogen states with the principal quantum number 
up to 15 and orbital quantum number up to 6 must be taken into account. The numerical 
procedure becomes unstable for 7 — » and 7—7-00. This non-convergency is caused by 
the fact that these points, where the symmetry changes (to 2D and ID, respectively), are 
peculiar for the perturbation theory. The dimension change causes the levels' degeneration, 
when a very large (divergent) number of levels is mixed due to perturbation, and has to be 
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FIG. 1: The maximum principle quantum number N of the states used in numerical calculations 
of the ground state exciton energy within a relative margin of 10 -4 as a function of the anisotropy 
parameter 7 (7 < 1). Solid line shows the logarithmic approximation for N. 

taken into account. To calculate the ground state exciton energy within a relative accuracy 
of 10 -4 , the levels with principle quantum number n < N should be taken into account. In 
Fig. 1 we show the numerically found dependence of N on the anisotropy parameter 7 (for 
7 < 1), which is approximately logarithmic. 

Note that at a rather strong anisotropy, when 7 -C 1, the ground state exciton behaves 
as E w -4 + IO.37 1 / 3 , |0 O (O)| 2 oc 7- 1 / 3 and yfjjF) oc 7- 1 / 3 , (see in Ref. EJ). Thus, it is 
useful to plot physical values in dependence on 7 1 / 3 instead of 7. 

A. Energy levels 

Figures 2 and 3 show the calculated eigenvalues X u of Eq. (|2"T|) . related to the exciton 
energies, E u = — 1/A 2 ,, as functions of 7 1 / 3 for 7 < 1 (left panels); 7 _1 ^ 2 A l/ are shown as 
functions of 7" 1 / 3 for 7 > 1 (right panels). The multiplier 7 -1 / 2 in the latter case makes the 
effective Rydberg finite when fi± — > 00. The binding energies of m = even parity states 
and m — 1 odd parity states are shown, respectively, in Figs. 2 and 3. Starting at 7 = 1 
from = = n, all the eigenvalues with the same m and parity do not intersect when 
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FIG. 2: Fock eigenvalues X u of m = even parity states as functions of the anisotropy parameter 
7 < 1 (left panel), and 7~ 1 / 2 A Jy as functions of 7 -1 / 3 , 7 > 1 (right panel). Solid curves never 
intersect each other due to a small anticrossing between the levels. The eigenvalues of purely 2D 
exciton (left panel) and ID exciton (right panel) are shown by semicircles. A linear approximation 
of the ground and first excited state eigenvalues is plotted by dotted and dashed lines, respectively. 

7 changes (multiple anticrossings occur due to the interaction between states) and approach 
the ground state eigenvalue of ID exciton 3 ^ 7~ 1 / 2 A l/ — > 7~ 1 / 2 Aq D — > (Figs. 2 and 3, right 
panels), when 7 — > 00. In the opposite case of 7 — > all shown eigenvalues approach the 
ground state eigenvalue Aq D = 1/2 of 2D exciton (m = 0, Fig. 2, left panel) or the first 
excited state eigenvalue A 2D = 3/2 (m = 1, Fig. 3, left panel). As it is clear from Fig. 2, the 
ground state eigenvalue dependence is almost linear over 7 1 / 3 for 7 < 1, and 
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En 



(30) 



(l+ 7 l/3)2- 

The ground state which lies much lower than the excited states almost does not interact 
with the latter. However, for the first excited state this interaction becomes much more 
significant, and its energy dependence upon 7 1 / 3 deviates from the linear one (cf. with 
dashed line in Fig. 2). The ratio of the energy separation between the ground state and 
the first excited state to the exciton binding energy is shown in Fig. 4. Starting from 3/4 for 
3D-isotropic exciton (Ei S — E 2 s) / E iS decreases monotonously with change of 7 and vanishes 
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FIG. 3: Fock eigenvalues X u of m = 1 odd parity states as functions of the anisotropy parameter 
7 < 1 (left panel), and 7~ 1 / 2 A l/ as functions of 7 _1//3 , 7 > 1 (right panel). 



TABLE I: Exciton binding energies of several lower states calculated by means of Brillouin-Wigner 

perturbation method, compared with that taken from Ref. 111. 

7 1 / 3 IS 2S 2P 2P± 3S a 3ZV 3P 3P± 
Ref. FL 0.8 1.233 0.3151 0.3663 0.2823 0.158 0.1375 0.1653 0.1272 
This work 1.2327 0.3151 0.3664 0.2823 0.1374 0.1577 0.1652 0.1272 
Ref.FL 0.4 2.01 0.695 0.933 0.3612 0.394 0.265 0.496 0.2100 
This work 2.011 0.6832 0.9381 0.3615 0.2835 0.4141 0.4959 0.2107 



°The levels classification used in the present work differs from that of Ref. 
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when 7 — > or 7 — > 00. Thus, this quantity can be considered as a measure of the anisotropy 
of a system. Note that within the fractional dimensional model (Ei$ — E 2 s)/Eis grows up 
as 1 — [(D — 1)/(D + l)] 2 (D is the dimensionality). Thus, in the anisotropic model the 
transition from 3D exciton to 2D or to ID exciton differs completely from that of a system, 
in which the carriers localization in one or two dimensions becomes stronger and using of 
the fractional dimensional model is justified. 

Results of our calculation for several low levels reproduce Faulkner's calculations 11 with 
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FIG. 4: The energy separation between the ground and first excited states in units of the exciton 
binding energy vs 7 1 / 3 . 



a good accuracy (see Table 1). As compared to Faulkner, we calculate a large number 
of excited states (up to 100 for each parity and m considered); we calculate the excitonic 
parameters in the region of 7 < 1 as well as 7 > 1, thus covering all possible values of 
the anisotropy parameter. Note the difference between Faulkner's and our designations of 
3S and 3Dq states.— When the states are split off due to perturbation, we always label the 
states with larger oscillator strengths at 7 ~ 1 as S-state, thus establishing an order reversed 
to that among the states with m/0, within our notations (see also discussions in Sec. IIII CI 
and Fig. 5). Thus, at 7 < 1 the 3D level lies lower than 3S, contrary to the classification 
by Faulkner -1 The same situation holds if we consider the higher excited states. 
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FIG. 5: Classification scheme of the energy levels of the anisotropic exciton with general quantum 
number n < 4 in accordance with the spherical approximation, Eq. (|35p . 7 = 0.5. 

B. Spherical approximation 

Even in case of small anisotropy |e| <C 1, the exciton states are linear combinations of 
hydrogen states with different /. However, for small e the admixture of such states becomes 
rather small, and the accounting only for the spherically symmetric part of the perturbation 
proves to be very useful for understanding the evolution of levels. It is important that within 
such a spherical approximation, the anisotropic exciton problem is exactly soluble. 

In this section we consider the approximate solution of the anisotropic exciton problem 
in a form ^(r) = R(r)Yi m (9, (p), thus taking into account only diagonal in I parts of the per- 
turbation, Eqs. (j24p . ()25p . and neglecting the perturbation matrix elements mixing different 
spherical harmonics. 

In order to neglect I 7^ I' matrix elements, let us replace in the Schrodinger equation, 
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Eq. (jSj), the operator p 2 z by the operator Q, defined as 

QY lm (6, <p) = p 2 Qi m Y lm (0, <p), (31) 
Qim = J cos 2 6\Y lm (6,<p)\ 2 d 3 <l (32) 

[see also Eq. ( |25|) ]. Then, after the substitution 

^irfo-' ^t£q-> (33) 

which, in fact, corresponds to a (I, m)-dependent mass renormalization, we arrive at a 
symmetrical (unperturbed) Schrodinger equation with the solution 

Mr) = (1 + eQ^r 3/2 0ll (—^-] , (34) 



l + eQ Zm 
n 2 (l + eQ im ) 

in units of Eq.(J2]) and with the use of dilatation of z. 

One can easily see from Eq. (|34|) that in this spherical approximation the perturbation 
compresses (for e < 0) or dilates (for e > 0) the scale of a given hydrogen wave function 
by the factor 1 + eQi m , which is different for different spherical harmonics. Note, that 
the hidden hydrogen-like symmetry is broken within this spherical approximation, and the 
binding energies now depend on I and m. However, the spectrum Eq. ( 135]) still has a 
hydrogen-like dependence on the principle quantum number n. 

In Fig. 5 we show schematically the energy levels of anisotropic exciton, calculated via 
Eq. ( 1351) for 7 = 0.5, n < 4 and all possible / and m. Equation ( )35l) provides a correct qual- 
itative description of the levels evolution and is in agreement with the result of calculations 
presented in Sec. IIII Al in the vicinity of 7 = 1 (see Fig. 6). 

The accounting for / 7^ /' matrix elements (in case of small e) yields correct quadratic in 
e terms in the energies. The rational form of matrix elements Eqs. ff26l) and (1271) allows us 
to sum up the standard perturbation theory series in the second order. The calculated in 
the second order exciton binding energies of several lower levels are given in Appendix A 
[see Eq. (IA17I) ]. their dependence on 7 is also illustrated in Fig. 6 (dashed lines). 

C. Oscillator strengths 

Within the envelope function approximation, the relative oscillator strengths of dipole- 
allowed transitions f u are proportional to 10^(0) | 2 (see, e.g., in Ref. Bearing in mind 
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FIG. 6: Lower energy levels as functions of 7 1 / 3 , calculated by means of the perturbation method 
(solid curves), within the spherical approximation (dotted lines), and in the 2nd order perturbation 
theory approximation (dashed curves). 



the expansion of Eq. (|29|) and the fact that for the unperturbed states (p u (0) 7^ only for 
I = m = 0, we get 



2 = Pu_ 



(36) 



Figure 7 shows the calculated numerically oscillator strengths of lower S, Dq and Go- 
like states as functions of the anisotropy parameter. It is seen in Fig. 7 that the oscillator 
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FIG. 7: The anisotropic exciton oscillator strengths of lower S, Dq and Go-like states as functions 
of the anisotropy parameter 7 1 / 3 , calculated numerically and within the spherical approximation 
(in units of the ground state oscillator strength at 7 = 1). Within the spherical approximation, 
the oscillator strengths of Dq and Go-like states are vanishing. 

strengths of all shown states do not vanish at 7 = 1. Originated from the degenerate states 
of isotropic 3D exciton, the perturbed states become fixed linear combinations of the former 
even when the perturbation tends to zero. It can be explained as follows. The perturbation 
of a symmetry lower than the original Hamiltonian implies the existence of strictly definite 
combinations of basis functions for degenerate states when 7 = 1, while the symmetry of the 
unperturbed Hamiltonian allows an arbitrary choice of these combinations. At 7 ~ 1 S-like 
state is optically more intensive than D -\ike state. The picture changes drastically with the 
increase of anisotropy. Near 7 1 / 3 = 0.8 the oscillator strength of 3D state overcomes that 
of 3S one. For 7 1 / 3 < 0.8 the intensity of the 35" state collapses due to the interaction with 
4Dq state and then revives after interaction with higher levels. Moreover, the anisotropy 
increase leads to substantial growth of the oscillator strengths of higher excited states, such 
as 4S and 4D > making them optically significant. Similar situation takes place if 7 > 1 
(when a transition from 3D to ID exciton occurs). Such a redistribution of the oscillator 
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FIG. 8: Energies and oscillator strength of excited optically active states vs 7 1 / 3 (7 < 1). The area 
of a circle is equal to the oscillator strength, normalized to the ground state oscillator strength, 
which is taken constant for all 7, see a single circle on the bottom (ground state) curve. 

strengths between different states is due to multiple unticrossings between energy levels 
interacting with each other. This effect can be clearly seen in Fig. 8, where the area of a 
circle placed on the energy curve is proportional to the oscillator strength of a given excited 
state, normalized to the ground state oscillator strength. 

Within the spherical approximation, as it follows from Eq. ( )34l) . 

/,= (l + £fV, (37) 

where = 1/n 3 are the oscillator strengths of the isotropic exciton (in units of f{g)- The 
oscillator strengths of IS", 2S and 3S states calculated according to Eq. (I57|) are displayed 
in Fig. 7 by dashed lines. Note that the oscillator strengths of 3S state, calculated numer- 
ically and within the spherical approximation, do not coincide at 7 = 1, as the spherical 
approximation does not reflect correctly the symmetry violation in the vicinity of this point. 
However, the sum of the oscillator strengths of 3S and 3D levels is equal to 1/3 3 . 
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IV. CONCLUSIONS 



The perturbation theory of anisotropic exciton is developed based of the Fock transfor- 
mation. This transformation depends on the exciton energies as adiabatic parameters and 
admits a separation of bound and scattering exciton states. For the bound states the eigen- 
functions are expanded into a complete set of hyperspherical harmonics on a 4D-sphere, 
creating a representation of the full symmetry group 0(4) of hydrogen-like system, and the 
perturbation matrix elements acquire an explicit algebraic form. This allows us to analyti- 
cally perform a partial diagonalization of the Hamiltonian matrix. It results in a spherical 
approximation which proves to be very useful for levels evolution analysis. The eigenvalues 
and eigenvectors are found by a numerical diagonalization of the effective Hamiltonian ma- 
trix. The energies and oscillator strengths of anisotropic exciton states are calculated for 
all values of the anisotropy parameter < 7 < 00 (including both flattened and elongated 
excitons), except the vicinities of 7 = and 7 = 00 where the dimensionality of the sys- 
tem changes, respectively, to D = 2 and to D = 1. It is found that with the increase of 
the anisotropy a strong redistribution of oscillator strengths between optically active and 
formerly inactive states occurs: the oscillations in optical intensities of higher excited states 
take place, and the switching on of formerly weak optical transitions is predicted. 
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Appendix A: Perturbation matrix. Exciton binding energies in the second order 
approximation 

Due to a separation of variables, the matrix element V ss >, Eq. (123]) . takes the form 




(Al) 
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where 



JJT' = Smm'MmM'm / if (x) P™ (x) X 2 dx 



-1 



Mi 



(2/ + 1) (Z-H)! 



1/2 



2 (Z + |m|)! 



7r(n + /)! 



Using the recurrent relations 



(A2) 

(A3) 
(A4) 

(A5) 



(2Z + l)a;if = (Z - m + + (Z + m)if 1? 



(A6) 



2(i/ + a)xC u Q (x) = (a + l)C£ +1 (z) + (2i/ + a - ljC^z), (A7) 
and normalization property for Legendre (if) and Gegenbauer (C^) polynomials we get 

1 - 4m 2 



rjm'i 

<~>IV 



1 + 



(2Z-l)(2Z + 3) 



Si v + 



1 

2Z + 3 



1 

21 - 1 



[Z 2 -m 2 ][(Z-l) 2 -m- 



(2Z + l)(2Z-3) 



H-2 v 



[(/ + l) 2 -m 2 ][(Z + 2) 2 -m 2 ] 
(2/ + 1)(2Z + 5) 



1+2 I' 



>(A£ 



and 



-a r 1 / (n-Q(n + I + l) 
X ""' " ^ + 2 V ^+1) ' 



1 / ( n + Z)(n-Z-l) , 

n+l n' T —\ 7 T7 On-1 n' 

2 y n(n — 1) 



To derive matrix elements Z^ 2 we use the tabulated integral 24 



7T 



2 r(i/ 



7r F(y) 



-C u m (l), m<n, 



where 



(m + 1v — 1 V 2 
(2z/ — l)!m! m 



^ 2 = P&ZWa £(1 " * 2 )(1 + ^(l - ^'-l^ttr^C^+i^)^- 
Using Eq. (1A7|) and the recurrent relations 



2i/(l - x 2 )C^(x) = (a + 2i/ - l)xC^_ x (x) - aC£(x) 
aCT^x) = 2(1/ - 1) K^x) - C£_ 2 (x)] , 



(A9) 



(A10) 

(All) 
(A12) 

(A13) 
(A14) 
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we are able to write the integral in Eq. (1A12|) in the form of Eq. (lAlOj) , 

?nn? = V; d V nll J^ j\l+x) l -*{l-xf-l (UC^-VmC 1 ^) (C^-Cl^dx, 

(A15) 

where 

= ( n + l){n + l-l) = (n-l)(n-l + l) 

Ann' " Ann' 

After simple transformations we arrive to Eqs. (I24l) -(l271). 

Note that we are able to apply the standard Reley-Schrddinger perturbation theory to 
Eq. ( 12 ip and to calculate analytically the perturbation theory corrections (to nondegener- 
ate levels) of a given order, due to the rational form of the perturbation matrix elements 
Eqs. ( I24p - (|27p . For example, for several lower levels the accounting for the perturbation 
theory terms up to the second order inclusive leads to 

E 1S = \ E, s = ! 

[i + ^+(i-ii)< 2 ] 2 ' 4[i+i e +(-^ + if)^] 2 ' 

EoPn — n, E<2p, = rT.(A17) 

4 h + A e+ _ 25441) £ 2l 2 ± 4 h + (8^ _ 49307 X 2 1 2 

L ~ 10 ~ V 35 21000/ J L ~ 10 ~ V105 63000/ J 
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